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ABSTRACT 

^ ■ We present new high and low resolution spectroscopic and photometric data of nine 

' members of the young association CMa Rl. All the stars have circumstellar dust at 

, some distance as could be expected from their association with reflection nebulosity, 

m ■ Four stars (HD 52721, HD 53367, LkHa 220 and LkHa 218) show Ha emission and 

we argue that they are Herbig Be stars with discs. Our photometric and spectroscopic 
, observations on these stars reveal new characteristics of their variability. We present 

first interpretations of the variability of HD 52721, HD 53367 and the two LkHa stars 
in terms of a partially eclipsing binary, a magnetic activity cycle and circumstellar 
■ dust variations, respectively. The remaining five stars show no clear indications of Ha 

Oh' emission in their spectra, although their spectral types and ages are comparable with 

those of HD 52721 and HD 53367. This indicates that the presence of a disc around a 
Jl^ ' star in CMa Rl may depend on the environment of the star. In particular we find that 

\ all Ha emission stars are located at or outside the arc-shaped border of the H ii region, 

^ . which suggests that the stars inside the arc have lost their discs through evaporation 

by UV photons from nearby O stars, or from the nearby (< 25 pc) supernova, about 
1 Myr ago. 



Key words: circumstellar matter - stars: early type - stars: evolution - stars: pre- 



' main sequence - stars: variables - open clusters and associations: CMa Rl 



1 INTRODUCTION 

Herbst & Assousa (1977) and Herbst, Racine & Warner 
(1978) found convincing evidence for a scenario in which 
the formation of the CMa Rl association was induced by 
a supernova explosion about 1 Myr ago. More recently, 
Comeron, Torra & Gomez (1998) argued that the star for- 
mation was triggered in pre-existing clouds and that it was 
only accelerated by the supernova shock. 

In a preceeding paper (Shevchenko et al. 1999; paper 
I) we have confirmed that most stars of CMa Rl have not 
been formed at the same occasion. Our objective of this 

* Based on observations collected at the European Southern Ob- 
servatory, La Silla, Chile 

t E-mail contact: Herman Tjin A Djie (herman@astro.uva.nl) 



paper is to present and discuss new optical and ultraviolet 
spectroscopic and photometric data for nine of its members. 
A tenth object, Z CMa (HD 53179), has also been observed, 
but since the interpretation of these data is more compli- 
cated, we will present the results separately (paper III). 
Apart from LkHa 220 (B5e) and LkHa 218 (B7e), which 
are close to ZAMS, the observed stars are of spectral type 
early B and on or above the Main Sequence (MS) with ages 
much larger than 1 Myr (paper I). It therefore seems useful 
to inspect their spectroscopic and photometric features in 
order to detect common characteristics in their atmospheric 
and circumstellar envelopes. 

While most of the stars have been observed before pho- 
tometrically and spectroscopically (e.g. Herbst et al. 1978; 
Finkenzeller & Jankovics 1984; Hamann & Persson 1992), 
the present optical spectra have a resolution which is seven 
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Figure 1. Observed (squares) and extinction-corrected (circles) spectral energy distributions of the nine programme stars. The solid 
lines show Kurucz (1991) models for the stellar photosphere, fitted to the extinction-corrected SED. 



times larger than in the earlier observations and for several 
stars a large amount of photometric observations has been 
made during the last decades. In addition we have surveyed 
information from the high and low resolution UV spectra, 
obtained by the International Ultraviolet Explorer (lUE). 
In section 2 we give a survey of the observations. Section 3 
describes the results of these and earlier observations. The 
photometric and spectral variability data of the four emis- 
sion line stars are discussed and interpreted in section 4. 
Section 5 summarizes the results and in this section we also 
make an attempt to understand these results in terms of the 
location of the stars in the CMa Rl association. 



2 OBSERVATIONS 

High resolution profiles of Ha, NalD, He I (5876 A), CallK 
and two lines of the Call (2) red triplet for HD 52721 and 
HD 53367 and of Ha, Na l D and He I (5876 A) for HD 52942, 
HD 53623, HD 53755, HD 53974 and HD 54439 have been 
secured in 1994, 1995 and 1996 at the European Southern 
Observatory (ESO), La Silla, Chile with the 1.4 m Coude 
Echelle Spectrograph (CES) on the Coude Auxiliary Tele- 
scope (CAT). Low resolution spectra of HD 52721, HD 53367 



and LkHa 220 were taken in 1992 and 1993 with the BoUer 
and Chivens spectrograph on the 1.5 m telescope at ESO. 
All spectra were reduced with the usual steps of bias sub- 
traction, fiatfielding, background subtraction and spectral 
extraction, and wavelength and flux calibration. Additional 
high resolution optical spectra of HD 52721 and HD 53367 
were obtained in October 1989 and November 1991 at the 
Special Astrophysical Observatory (SAO) of the Russian 
Academy of Sciences with the Main Stellar Spectrograph 
(MSS) on the 6-m reflector telescope. These spectra were 
processed on the automatic spidometer and microdensito- 
meter at the Crimean Astrophysical Observatory. Details of 
all observations are given in Table 1. 

Photometric observations of the stars have been made 
at ESO in various pass-bands {UBVRIJHKLM and Wal- 
raven WULBV; optical photometry not published before 
is given in Table 5) which, together with the IRAS point 
source fluxes (paper I) and UV fluxes from the ANS and 
TDl catalogues (Wesselius et al. 1982; Thompson et al. 
1978), permitted us to determine the spectral energy dis- 
tribution (SED) of each star. For HD 52721 (GU CMa), 
HD 53367 (V570 Mon), LkHa 218 (HT CMa) and LkHa 220 
(HU CMa) we also constructed the light curves and colour- 
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Table 1. Log of spectroscopic observations. 



ESO 



Star 


Rang 


;e [A] 


Date 


JD+2400000 


R [mA] t 


cxp [min] 


HD 52721 


4200- 


-7800 


11 Dec. 92 


48967.7653 


1076 


2 




6522- 


-6602 


11 Dec. 96 


50428.2279 


65.6 


15 




5849- 


-5920 


13 Dec. 96 


50430.1797 


58.8 


15 




3912- 


-3955 


14 Dec. 96 


50431.3548 


39.3 


15 




8466- 


-8572 


15 Dec. 96 


50432.2681 


85.2 


20 


HD 53367 


4200- 


-7800 


11 Dec. 92 


48976.7722 


1076 


2 




6536- 


-6590 


16 Dec. 94 


49702.3472 


65.6 


15 




6540- 


-6590 


15 Jan. 95 


49732.2750 


65.6 


20 




5860- 


-5910 


16 Dec. 94 


49702.2125 


58.8 


20 




5860- 


-5910 


16 Jan. 95 


49733.2887 


58.8 


15 




3912- 


-3955 


14 Dec. 96 


50431.3561 


39.3 


15 




8466- 


-8572 


15 Dec. 96 


50432.2354 


85.2 


20 


LkHa 220 


6400- 


-8950 


1 Jan. 94 


49353.2103 


1245 


25 




3600- 


-6800 


2 Jan. 94 


49354.3145 


1563 


15 


HD 53623 


5849- 


-5920 


13 Dec. 96 


50430.1921 


58.8 


20 




6522- 


-6602 


12 Dec 96 


^0429 2371 


65.6 


20 


HD 53755 


5849- 


-5920 


13 Dec. 96 


50430.1921 


58.8 


20 




6522- 


-6602 


12 Dec. 96 


50429.2013 


65.6 


15 


HD 52942 


5849- 


-5920 


14 Dec. 96 


50431.1600 


58.8 


15 




6522- 


-6602 


12 Dec. 96 


50429.2208 


65.6 


20 


HD 53974 


5849- 


-5920 


14 Dec. 96 


50431.1445 


58.8 


20 




6522- 


-6602 


12 Dec. 96 


50429.2139 


65.6 


7 


HD 54439 


5849- 


-5920 


14 Dec. 96 


50431.1719 


58.8 


10 




6522- 


-6602 


12 Dec. 96 


50429.1852 


65.6 


20 




5849- 


-5920 


14 Dec. 96 


50431.1807 


58.8 


25 


SAO 


Star 


Rang 


;e [A] 


Date 


JD+2400000 


R [mA] 




HD 52721 


3700- 


-4900 


14 Oct. 89 


47814. 


200 






3700- 


-4900 


15 Oct. 89 


47815. 


200 




HD 53367 


3700- 


-4900 


15 Oct. 89 


47815. 


200 






3700- 


-4900 


02 Mar. 91 


48318. 


200 




lUE high res. 


Star 




SWP 


Date 


LW 


Date 



HD 52721 10592 12 Nov. 1980 R09283 12 Nov. 1980 

HD 53367 38686 27 Apr. 1990 P17816 27 Apr. 1990 

HD 53755 27974 21 Mar. 1986 R02166 24 Aug. 1978 

HD 53974 10580 10 Nov. 1980 R09272 10 Nov. 1980 



magnitude diagrams from the photometric sequences ob- 
tained by the observers at Van Vleck Observatory (Herbst 
et aL 1994), at Mt. Maidanak Observatory (Kilyakov & 
Shevchenko 1976) and at Corralitos Observatory (Halbe- 
del 1989, 1991, 1999). High and low resolution ultraviolet 
spectra were taken from the Final Reduction Archive of the 
WE (processed by NEWSIPS) at the ESA Satellite ground 
station in Villafranca del Castillo (Madrid). 



3 SPECTROSCOPIC RESULTS 
3.1 The Photospheres 

Spectroscopic classifications from the visual spectral range 
of the nine stars have been published frequently and are col- 
lected in Table 2. For the early B-type stars the photospheric 
information from the visual spectrum is limited to some lines 



Table 2. Properties of programme stars. 



Star 


Var. Name Sp. T. 




Vj 


V sin i 










[pc] 


[mag] 


[km s"-'] 


[km s"-'] 


HD 52721 


GU CMa 


Bl Ve 


600 


6.40-6.80 


350-450 


-f9-+22 






BllVe 


800 








HD 53367 


V750 Mon 


BO IVe 


775 


6.90-7.20 


25-50 


-1-11 — 1-24 






BO Ille 


1050 








LkHa 218 


HT CMa 


B6Ve 


900 


11.61-12.05 


50 




LkHa 220 


HU CMa 


B5Ve 


1090 


11.60-12.19 


75 




HD 52942 


FZ CMa 


B2IV 


950 


8.05-8.44 


370 




HD 53623 




Bl V 


1150 


7.93-8.01 




-1-30 


HD 53755 


V569 Mon 


B0.5IV 


1010 


6.40-6.53 


410 




HD 53974 


FN CMa 


BOH 


1025 


5.38-5.43 


150 


-H31 


HD 54439 




Bl V 


960 


7.68 


260 





of Mgll and He I, which may have resulted in some uncer- 
tainties in the spectral classifications. In the UV, however, 
more purely photospheric lines can be found: Hell (1640 A), 
Cm (1247 A), Niii (1747, 1751 A), Sim (1299, 1417 A) and 
many lines of Felll, Feiv and Fe v (see e.g. Prinja 1990). In 
the lUE archive high resolution spectra of four members of 
CMa Rl are available: HD 52721, HD 53367, HD 53755 and 
HD 53974. We have compared these spectra with those in 
the UV spectral atlasses for B star classification by Rountree 
& Sonneborn (1993) and Snow et al. (1994). 

From this comparison we derive the types Bl IV-V for 
HD 52721 (comp. stars: HD 37303 BIV with usini = 260 
km s"^ and HD 108483 B2V with vsini = 245 km s"^), 
B0IV,V for HD 53367 (comp. stars: HD 36512 BOV with 
vAni = 20 km s"^ and HD 34816 BO. 5 IV with vsmi = 
67 km s^^) and B0.5III,IV for HD 53755 and HD 53974 
(comp. stars: HD 219188 BO. 5 HI with vsmi = 265 km s"^ 
and 5 Sco (HD 143275) B0.3IV with «sini = 180 km s"^). 
In this procedure approximate values of v sin i = 400, 30, 410 
and 150 km s~^ have been taken into account for the stars 
HD 52721, HD 53367, HD 53755 and HD 53974 respectively 
(Uesugi & Fukuda 1976; Finkenzeller 1985; Bohm & Catala 
1995). From the high ionization lines in the comparison spec- 
tra it also becomes clear that HD 53755 and HD 53974 have 
fast outflows. Our classifications from the UV spectra are 
close to those of the spectral and luminosity types from the 
visual spectra which were obtained earlier. Note that when 
we use the My values of Schmidt-Kaler (1982), these classi- 
fications are consistent with a distance of --^ 1 kpc (close to 
the value of 1050 pc, derived for CMa Rl in paper I). 

3.2 The Outer Circumstellar Regions 

Let us now inspect the information on the circumstellar dust 
and gas of the nine stars. Fig. 1 shows the spectral energy 
distributions (SEDs) of the stars. These distributions are fit- 
ted by the calculated Kurucz (1991) model distributions for 
the spectral and luminosity class of each star, together with 
the observed E{B — V). In order to get a good agreement in 
the visual and UV interval, one has to vary the extinction 
law, characterised by the value of Rv (the ratio of total ab- 
sorption in the photometric V-band over the colour excess 
E{B — V)). The best fitting values of Rv, as well as the 
stellar parameters used for these fits, are included in Fig. 1. 
The colour excess of HD 53367 is exceptionally large 
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Figure 2. Low resolution /(7J? spectra of 6 stars. Tlie arrows near 
the stellar names point to the corresponding flux scales. 



compared to those of the other stars in CMa Rl. This is 
also expressed by the 2200 A absorption in the low resolu- 
tion UV spectrum of this star, which is much stronger then 
those in the corresponding spectra of the other stars (Fig. 2). 
From Fig. 1 it is clear that only for the four Ha emission 
line stars HD 52721, HD 53367, LkHa 218 and LkHa 220 the 
extinction corrected SED has an excess flux at wavelengths 
longward 1.58 fim. Longward about 4 ^m, the excess for 
LkHof 218 and LkHa 220 increases strongly and can only 
be explained by thermal re-emission from nearby hot cir- 
cumstellar dust. For HD 52721 and HD 53367 the infrared 
excesses longward 4 nm are smaller and could be due to the 
presence of cooler circumstellar dust at larger distances from 
these stars, or to the contribution of cool stellar companions 
(see Sect. 4.2). The excesses at shorter wavelengths are much 
smaller and could be contributed by free-free emission in the 
inner circumstellar regions (see Sect. 3.3). 

In order to find more information about the gas com- 
ponent of the circumstellar medium of the stars we will first 
discuss the high resolution profiles of the NalD (Fig. 3) 
and CallK lines (Fig. 4c). For HD 52721 we find equal dis- 
placements (4-19 km s"^) and FWHM's (22 km s~^) for 
the CallK and NalD lines. For HD 53367 we find similar 
values (displaced -|-17 km s~^ and 24 km s~^) for CallK 
and NalD. This indicates that the resonance line of Can 
is formed in the same region as those of Nai. The shape of 
the Na I D profiles (Fig. 3) shows the presence of more than 
one component, which could be contributed by interstellar 
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1 

ll 







5885 



5890 5895 
A [A] 



5900 



Figure 3. NalD profiles of 7 stars. The horizontal and vertical 
dashed lines show the estimated continuum level and the Na I D2,i 
rest wavelengths, respectively. 



clouds in the CMa association or else from circumstellar gas. 
The redshifts of HD 52721 and HD 53367 are within the 
range of their radial velocity variations (Abt & Biggs 1972). 
The main components are redshifted by about 19 km s~^ 
for HD 53974, HD 53755 and HD 54439, which is within 
the range of the radial velocities measured for the double 
stars HD 53974 and HD 53755 (Abt & Biggs 1972). For 
HD 53623 we measure -1-12 km s~^ for the main (blue) com- 
ponents and +27 km s~^ for the second component. The 
radial velocity from the literature (Abt & Biggs 1972) is 
+30 km s"\ There are no indications that this star is a bi- 
nary. HD 52942 is a double line eclipsing binary with equal 
mass companions (Moffat & Vogt 1974), but we see only 
one NalD doublet in our (noisy) spectrum. The difference 
in redshift may be hidden by the noise, but the average value 
(+15 km s~ ) can only be understood if the absorption is 
formed at some distance of the binary system. From the 
general agreement between the NalD redshifts and the stel- 
lar radial velocities we conclude that the NalD doublet is 
formed in the (wide) circumstellar region of the stars. From 
the EWs of the NalD doublet lines we can evaluate the 
column density of the foreground gas with the Stromgren 
(1948) doublet-ratio method. In Table 3 we give the EWs of 
the NalD lines and the derived column densities of NalD. 
From A'^(Nal) and by assuming a galactic abundance ratio 
7V(H)/A'^(Nal) of 5 x 10^ we can derive colour excesses for 
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Table 3. Comparison of interstellar colour excess, derived from the Nal{l) doublet with observed total colour excess. 



Star 


EW (NaDi) 


EW2/EW1 


logAf(Na) 




E(B - V)ohs 




u.ouo 


1 32 


12 72 


n 1 zL^ 
u. ±^o 


n 9S 

U.Zo 


HD 52942 


0.393 


1.18 


13.05 


0.222 


0.42 


HD 53367 


0.461 


1.17 


13.18 


0.259 


0.72 


HD 53623 


0.375 


1.30 


12.84 


0.170 


0.22 


HD 53755 


0.311 


1.24 


12.85 


0.171 


0.21 


HD 53974 


0.410 


1.22 


13.07 


0.226 


0.34 


HD 54439 


0.399 


1.21 


13.00 


0.207 


0.31 


LkH« 218t 


0.117 


1.23 


12.42 


0.098 


0.57 


LkHa 220t 


0.141 


1.03 


13.71 


0.513 


0.43 



^The EWs for LkHa 220 and LkHa 218 are less accurate than those for the other stars, because they were measured from the figures in 
the paper of Finkenzeller & Mundt (1984). 



an average interstellar dust/gas ratio with the empirical re- 
lation of Hobbs (1974): Ar(Nal) = 1.7 X lO^"* x EiB-Vf "^. 
The results are listed in column 5 of Table 3 and can be com- 
pared with the values of E{B — V) (column 6), derived from 
the photometric observations by fitting the SEDs. From this 
comparison it is clear that the interstellar dust gives signif- 
icantly lower colour excesses than those obtained directly 
from the photometry (see also Paper I Sect. 3.1). This can 
be understood if the dust/gas ratio in the circumstellar and 
intra-association cloud regions is higher than in the general 
interstellar foreground. From Table 3 we see that this ratio 
will be especially high for HD 53367. 

We have inspected the high resolution lUE spectra of 
HD 52721, HD 53367, HD 53755 and HD 53974. In these 
spectra one finds a large number of narrow lines with zero ex- 
citation energy of their lower levels, which will be formed in 
the interstellar foreground; e.g. at 2851 A (Mgl (1)), 2026 A 
(Mgl(2)), 2795 A, 2802 A ((Mgll(l)), 2593 A (Mnll(l)), 
2585 A (Feii(l)), 2382 A (Feii(2)), 1526 A (Sin (2)). In 
contrast, the number of low excitation lines in these spectra 
is small and these lines are very weak. We note that only 
a few lines (2550.8 A of Fen (240), 2556.6 A of Mn(20), 
2666.6 A of Fell (263) and 2724.5 A of Mnii (33)) with lower 
levels up to ^ 3.5 eV are present in the spectra of HD 52721 
and HD 53974 but are missing in those of HD 53367 and 
HD 53755. Since the latter two stars are not far apart from 
each other in CMa Rl, their common lack of these low exci- 
tation lines could be due to the conditions in their common 
foreground interstellar clouds. The scarceness and weakness 
of the low excitation lines indicate that the four stars which 
we inspected in the UV have only a very small amount of 
circumstellar gas, which is consistent with the high dust /gas 
ratio in the circumstellar region. Only one of the stars, 
HD 53367, shows the possible presence of the Mgll (3) dou- 
blet (2790.77 and 2797.99 A) in one of its spectra. 

3.3 The Inner Circumstellar Regions 

3.3.1 The B-emission line stars: HD 52721, HD 53367, 
LkHa 218 and LkHa 220 

Of special interest are the emission lines of Ha (Fig. 4a), 
the red Ca ll (2) lines (Fig. 4d) and the emission compo- 
nent in the He I lines in the visual range of the spectra of 
HD 52721 and HD 53367 (Fig. 4b). In high resolution. Ha 
shows emission profiles (Fig. 4a) which are symmetrical for 



HD 52721 and asymmetrical (distorted by a red-shifted ab- 
sorption component) for HD 53367. For the first star the 
shape does not vary much on short (days) and long (yrs) 
timescales (see Sect. 4). For HD 53367 small shape vari- 
ations in the red wing occur on time-scales of one month 
or less (Fig. 4a) and large shape variations over the whole 
profile occur over time-scales of years. Because of its very 
low vs'mi {~ 30 km s~^) we could correct the Ha profile of 
HD 53367 for its underlying photospheric profile (assumed 
to be that of a Kurucz model Ha profile for T^b = 30,000 K 
and \ogg = 3.5) and conclude that the outflow velocity is 
100-150 km s-^ 

For HD 52721 with its large wsini (350-400 km s"^) 
we have no comparison profile, but the observed profile in- 
dicates that the outflow velocity of this star is less than 100 
km s~^. The parameters of Ha of the four stars are given 
in Table 4, which shows that they are variable in time. For 
HD 53367 we measured an EW of —6.1 A, which strongly 
differs from -14.0 A (Andrillat 1983) and -16.2 A (Hamann 
& Persson 1992), but is close to the —5.0 A, measured by 
Corcoran & Ray (1998) for their proflle of Jan. 1, 1992. We 
measure a FWHM of 260 km s"^ for HD 52721 and 280 
km s~^ (or 310 km s~^ if we account for the presence of an 
absorption component) for HD 53367. For LkHa 218 and 
LkHa 220, emission proflles of Ha on Nov. 15, 1981 have 
been presented by Finkenzeller & Mundt (1984). These pro- 
flles have double peaked shapes with blue-shifted (« —70 
km s"'^) absorptions. The variability of Ha is further dis- 
cussed in Sect. 4. 

For HD 52721 and HD 53367 we have derived volume 
emission measures from the flU-in of the Balmer discontinu- 
ity by bound-free and free-free emission in the circumstel- 
lar envelope. This flU-in could be determined by comparing 
the Balmer jump, derived from photometric observations on 
the Walraven system (which has two pass bands shortward 
and two passbands longward the Balmer limit at 3704 A), 
with the model-predicted Balmer discontinuity (e.g. from 
the non-LTE models of Mihalas). The observed Walraven 
magnitudes and the derived Balmer flU-in ADf, are given 
in Table 5. With the method of Garrison (1978) and the 
assumption that Tea ~ 10,000 K, we then derive volume 
emission measures E = J neUidV/R^ of 6.1 x 10''^ cm~^ for 
HD 52721 and 1.3 x 10^^ cm"'' for HD 53367. If we assume 
an electron density of 10^^ cm"'' in the circumstellar enve- 
lope region close to the stars we estimate with the formula 
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Table 4. Equivalent Widths (EWs) and Full Width Half Max. (FWHMs) of four emission-line stars in CMa Rl. 



Star 


A [A] 


Line 


EW [A] 


Error [A] 


FWHM [k 


m s ^] Date 


Ref. 


HD 52721 


3933.6 


CallK 


+0.14 


0.03 


25 


14 Dec. 


1996 


This paper 




5875.6 


He I 


+0.67 


0.10 


490 


13 Dec. 


1996 


This paper 




6562.8 


Ha 


-5.2 






330 


14 Nov. 


1981 


Finkenzeller & Mundt (1984) 




6562.8 


Ha 


-9.4 








31 Jan. 


1983 


Andrillat (1983) 




6562.8 


Hq 


-10-14 






300 


4— y Mar. 


1985 


Praderie et al. (1991) 




6562.8 


Ha 


-8.8 








23 Dec. 


1991 


Corcoran & Ray (1998) 




6562.8 


Ha 


-14 








11 Jan. 


1995 


Oudmaijer & Drew (1999) 




6562.8 


Ha 


-9.9 


0.15 




11 Dec. 


1996 


This paper 




7773 


Ol 


-0.24 






311 


Dec. 


1987 


Hamann & Persson (1992) 




8446.5 


O I 


—2.06 








Dec. 


1987 


Hamann & Persson (1992) 




8465 


P17 


-0.35 






335 


22 Jan. 


1983 


Andrillat et al. (1988) 




8498.02 
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Figure 4. High resolution profiles for HD 52721 and HD 53367. (a) Ha, (b) Hel (5876 A), (c) CallK and (d) Call(2). For HD 53367 
the black and grey profiles of Ha and Hel correspond to the two exposures, one month apart (see Table 1). 



Table 5. Walraven photometry of HD 52721 and HD 53367 and 
corresponding fill-in of Balmer decrement, emission measure and 
mass loss rate. 





HD 52721 


HD 53367 


Date 


May 1990 


25 Jan. 1981 


Vw 


0.153 ± 0.002 


-0.002 ± 0.002 


(V-B)w 


0.036 ± 0.002 


0.202 ± 0.002 




0.018 ± 0.002 


0.055 ± 0.002 


{B-U)w 


0.016 ± 0.003 


0.043 ± 0.003 


{B-W)w 


0.027 ± 0.003 


0.146 ± 0.003 


Vj [mag] 


6.48 


6.90 


ADfc [mag] 


0.24 


0.30 


EM [cm^l 


6.2 X lO^*^ 


1.3 X 10^^ 


VHa [km s-1] 


Ri 100 


140 


M [Mo yr-i] 


9.2 X 10"^ 


4.2 X 10"'' 



of Garrison (1978) that the mass outflow in the circumstel- 
lar envelope close to the star is 9.2 x 10"'' Mq yr~^ for 
HD 52721 and 4.2 x 10"® M© yr"^ for HD 53367. 

The high resolution profiles of the emission lines at 
8500 and 8543 A (Fig. 4d) allow us to separate the 8498.02 
and 8542.09 A lines of the Call (2) triplet from the partly 



blended Paschen lines P16 (8502.49 A) and P15 (8548.39 A). 
As a result of the blending, the peaks are slightly shifted. 
For HD 52721 we have decomposed the profiles into two 
components with a FWHM of 225 km s~^ each, which 
gives us EW(P16) = -0.34 A, EW(Caii 8498) = -0.46 A, 
EW(P15)= -0.37 A and EW(Call 8542) = -0.75 A. To- 
gether we have -0.80 A for the 8500 profile and -1.12 A for 
the 8543 profile, compared with —0.72 A and —0.55 A de- 
termined by Andrillat, Jaschek & Jaschek (1988). Note that 
these authors give an EW of —0.35 A to both unblended 
P14 (8598 A) and P17 (8465 A) lines. We therefore may ex- 
pect a similar EW for P15 and P16, which is very close to 
our results above. The EWs of the blended lines quoted by 
Andrillat et al. differ from our values in Table 4, which may 
indicate a variability of the Call (2) fluxes. 

For HD 53367 we made similar decompositions with 
FWHM values of 205 km s'^ This gives EW(P16) = 
-0.56 A, EW(Caii 8498) = -0.83 A, EW(P15) = -0.74 A 
and EW(Call 8542) = -1.16 A. Andrillat et al. (1988) give 
EW(P14) = -1.54 A and EW(P17) = -1.11 A for this star. 
This is consistent with our result above that the EW of Ha 
in 1992 was very low (much smaller than in Oct. 1981). The 
blended lines at 8500 and 8543 A were observed by Andrillat 
et al. in Jan. 1983. Our most important conclusion from the 
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Figure 5. Low resolution spectra of LkHo 220, HD 52721 and HD 53367. Insets sliow Ha and the red spectrum of LkHa 220 in more 
detail. Lines, dots and stars in the inset mark the positions of lines of Hi, Call and Ol, respectively (see Table 4). 



high resolution profiles of HD 52721 and HD 53367 is that 
the emission of the Call (2) lines is significant and has a 
FWHM of the order of 210 km s~\ Parameters of the Call 
lines are collected in Table 4. 

Since Call cannot exist at the stellar photosphere the 
emission of Call (2) must originate from the envelope region. 
Hamann & Persson (1992) have given arguments for the 
presence of discs around Herbig Ae/Be stars with Ha and 
Can emission. For early B-type stars these discs are proba- 
bly very similar to those of classical Be stars. The FWHM's 
of the Ca ll (2) lines suggest that Ca ll (2) is emitted in an 
outer region of a circumstellar disc (at a Kepler radius of 
29 R* for HD 53367 and 17.2 for HD 52721), whereas 
the Ha emission could come from a region of the disc which 
is somewhat closer to the star (at a Kepler radius of 13.5 R* 
for HD 53367 and 12.4 R* for HD 52721). Here we have 
implicitly assumed that both star-discs are being observed 
more or less edge-on. It is interesting to note that the low 
resolution red spectrum of LkHa 220 (Fig. 5) also shows the 
presence of strong emission in the Call (2) triplet. 

The profiles of the O I emission lines for HD 52721 and 
HD 53367 in the red part of the spectrum have been ob- 
served by Hamann & Persson (1992). The emission profiles 
of these lines are weak and broad (Fig. 3 in their paper) 
and only the widths of the blended 8446 A Ol (4) line could 
be estimated (FWHM of ~ 300 km s"^). The authors re- 
port EWs of -0.24 A and -0.37 A for the unblended triplet 
Ol(l) line (7773 A) of HD 52721 and HD 53367 respec- 
tively and -2.06 A and -2.68 A for the blended Ol(4) 
line (8446 A) of HD 52721 and HD 53367 respectively. The 
8446 A O I (4) emission is much stronger than that of the 



7773 A line, which cannot completely be explained by the 
blend with P18 (8437 A). This Ol(4) emission is common 
in classical Be stars and probably due to pumping of the 
^Do state of Ol by Ly/3 (Bowen 1947). The fact that the 
1302 A O I (2) line is in absorption, while for the same stars 
the O I (4) transition, which feeds the 1302 A transition, is in 
emission has led Oegerle, Peters & Polidan (1983) to argue 
that because the optical depth of line scattering in the res- 
onance 1302 A line is much larger than in the subordinate 
8446 A line, the 1304 A photons could be scattered out the 
line-of-sight in equatorially viewed flattened Be envelopes. 
It would suggest that the O I emission is coming from flat- 
tened discs around HD 52721 and HD 53367 and that both 
discs are viewed edge-on. This suggestion seems to be con- 
firmed by recent Ha polarization measurements (Oudmaijer 
& Drew 1999) of HD 53367 which indicate that this star is 
not observed pole-on. In the low resolution red spectrum of 
LkHa 220 (Fig. 5) we can distinguish the emission in O I (1) 
at 7773 A, and that in Ol (4) at 8446.5 A, slightly blended 
with P18. 

The He: (5876 A) absorption profile of HD 52721 (Fig. 
4b) is very broad with a FWHM of ~ 450 km s~^, approxi- 
mately the same as for the He I lines at 4026 and 4471 A in 
the spectrum of this star, obtained by Finkenzeller & Mundt 
(1984). This large width probably reflects the photospheric 
origin of these He I lines in this fast rotating star. Inside the 
line profile of the 5876 A line a weak emission component can 
be seen at both sides of the line centre. The He: 6678 A pro- 
files of this star (Corporon & Lagrange 1999) also have such 
an emission component, which shows variability in strength 
on a time-scale of less than one day. Such emission compo- 
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Figure 6. High resolution spectral intervals of HD 52721 taken 
on Oct. 14, 1989 with the Main Stellar Spectrograph on the 6-m 
BTA telescope of the SAO. The upper panel shows profiles of If7 
and He I (4387.9 and 4471.5 A). The lower panel depicts those of 
He, CallH and K and Hel (4009.3 and 4026.2 A). 



nents have been observed more often in the Hel hnes of B 
and Be stars and have been interpreted as CQE's (central 
quasi-emission peaks), reduced absorption which can occur 
when the steUar disc is occulted by circumstellar gas mov- 
ing in Kepler rotation around the star and when we observe 
the circumstellar disc edge-on (Hanuschik 1995). It has been 
emphasized by Rivinius, Stefl & Baade (1999) that such fea- 
tures develop in the inner regions of circumstellar discs and 
are an indication of gas transfer from the star to the disc. In 
HD 52721 the high rotational line broadening will prevent 
that such CQE's are observed in the weaker metal lines. For 
this star, however, an interpretation of the feature in terms 
of a tidal effect due to the nearby companion (see Sect. 4) 
should also be considered. This should be checked by taking 
spectra at different phases of the binary motion. 

The Hel lines of HD 53367 have photospheric absorp- 
tion profiles with an average FWHM of 150 km s~^. In the 
5876 A hue profiles (Fig. 4b) we observe the presence of 
an emission component, with a blue-shift of ~ 150 km s^^ 
with respect to the line centre. The profile of Jan. 16, 1995 
shows somewhat larger emission and deeper absorption com- 
ponents than that, taken 31 days before. These components 
are also clearly present in the corresponding profiles for 
HD 53367 and perhaps LkHa 220 obtained by Finkenzeller 
&: Mundt (1984) in Nov. 1981. The emission and absorption 
components of this line seem to vary in time (see Table 4). 
Inspection of the He I 6678 A profiles of HD 53367 (Cor- 
poron & Lagrange 1999) show variability on time-scales of 
less than 26 days, which could partly be due to radial veloc- 
ity variations with a probable (binary) period of 166 days 
but changes in the width of this absorption profile suggest 
also here the presence of a variable emission component. 
Because the binary parameters suggest a large separation 



of the companions the variations can not be due to tidal 
interactions. Since the emission feature in the profile of the 
5876 A line is not in the centre of the line we think it cannot 
be interpreted as a CQE, but the blue-shift of the emission 
suggests a hot plasma flowing from the photosphere to the 
inner circumstellar disc. We should add here that the profile 
of the photospheric He ll 1640 A line also shows indications 
for the presence of blue-shifted (« 150 km s~^) emission. 
The parameters of the proflles of the Ha, Hel (5876 A), 
NalD and Can (2) lines of HD 52721 and HD 53367 are 
collected in Table 4. 

The proflles of several other lines in the spectrum of 
HD 52721 (H7, Hei (3888.6, 3964.7, 4437.9, 4471.5 A), He, 
and CallH and CallK) are shown in Fig. 6. For H7 we find 
an EW of 3.60 ± 0.06 A, or 3.96 ± 0.04 A, if we correct 
for a small contribution of emission in the line centre. This 
can be compared with an EW of 3.0 A, estimated from the 
spectrum of Finkenzeller & Jankovics obtained on Feb. 18, 
1981. The value of the EW of H7 is consistent with the 
prediction of the Kurucz model for T^ff = 25,000 K and log g 
— 4.0 (Bl IV), which supports the classification of Table 2. 

In Fig. 5 we show the low resolution spectra of the Ha 
emission stars HD 52721, HD 53367 and LkHa 220. The lat- 
ter star has the strongest emission in Ha and H/3, but a good 
comparison is difficult, since the three stars are variable in 
emission. The blue low resolution spectra of HD 52721 and 
HD 53367 are similar to the spectra of these stars which were 
obtained by Finkenzeller & Jankovics (1984) with higher 
resolution. Their main absorption lines are those of Hel as 
expected for these early B spectral types. Comparison of the 
higher resolution blue spectrum of HD 53367 published by 
Finkenzeller & Jankovics (1984) with that of the spectrum 
of r Sco (BOV) in the same paper, shows faint absorption 
lines of On (2), (5), (6), (10) and a faint line of Siiv (1) 
at 4116 A, which also occur in the spectrum of r Sco, and 
therefore probably are formed in the photosphere. The pres- 
ence of such faint lines cannot be verified in the spectrum 
of HD 52721, because of the broadening of these lines in 
this fast-rotating star. In contrast to the blue spectra of 
HD 52721 and HD 53367, that of LkHa 220 (B5V) shows 
emission in the photospheric Fen (42) lines at 4923.9, 5018.4 
and 5169.0 A (BO with 30,000 K is too hot for Fe 11). The He i 
lines at 5876 and 6678 A and the NalD lines are dominated 
by absorption, as well as the Cal lines at 4026 and 4435 A 
and those of probably Mgll (4481 A) and Till (4469 A). 
There is a weak indication of emission of [Ol] (6300 A), 
also observed by Bohm & Catala (1995). The low resolution 
blue spectra of LkHa 220 and LkHa 218, taken by Cohen 
& Kuhi (1979) on Jan. 14, 1977 show no clear lines except 
the Balmer and NalD lines. In these spectra the Ha emis- 
sion of LkHa 220 is stronger than that of LkHa 218. In the 
red spectrum of LkHa 220 (Fig. 5) one can distinguish the 
7774 A and 8446 A lines of Ol(l) in emission. The EW 
of the unblended 7774 A line is ~ —0.6 A, compared with 
+0.6 A for a standard B5IV-V star (Fig. 2 of Slettebak 
1986). As for HD 52721 and HD 53367, the red triplet of 
Ca 11 (2) of LkHa 220 is in emission and blended with emis- 
sion from P13, P15 and P16 (see insert of Fig. 5). In this 
spectrum the emission fiux of the unblended Paschen lines 
is decreasing from a strong Pll emission at 8863 A to a very 
weak emission in the P17 line. 

In the high resolution lUE spectra most of the UV lines 
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Figure 7. High resolution spectral intervals of HD 53367 (spectra of the Main Spectrograph of the SAO). The left panels show the 
profiles of H/3, whereas the right panels display the profiles of H(5, He I (4120.8 and 4471.5 A) and Mgll 4481 A. The panels correspond 
with observation dates 02 March 1991, at maximum brightness (upper panels) and 15 Oct. 1989, at minimum brightness phase of the 
star (lower panels). 



are formed in the photospheres and the cooler parts of the 
inner circumstellar region. The profiles of the Siiv reso- 
nance lines have blue wings extending up to 530 km s~^ 
for HD 52721 and up to 470 km s'^ for HD 53367 (Fig. 8). 
Those of the Civ resonance lines extend up to 340 km s~^ 
for HD 52721 and to 370 km s"^ for HD 53367. These lines 
are formed in outflowing winds, accelerated in the circum- 
stellar envelope. 



3.3.2 The B stars: HD 52942, HD 53623, HD 53755, 
HD 53974 and HD 54439 

From the Walraven photometry (The, Wesselius & Janssen 
1986) we measured the flU-in of the Balmer discontinuity 
of HD 53755 and HD 53974 and found values of ADb of 
0.15 and 0.20 mag respectively. The lack of Hq emission in 
the former two stars must therefore be due to a much lower 
Tsheii (< 10,000 K) and consequently a much lower emission 
measure for the inner circumstellar envelope of the stars. 
This temperature may be too low to populate the upper 
level of the Balmer series. 

For the five B stars, in which Ha is in absorption, there 
is a certain similarity in profile shape between Ha (Fig. 9a) 
and He I (5867 A) (Fig. 9b). This suggests that Ha, similar 
to He I (5867 A), is formed in the photosphere of these stars. 
The profile parameters (EW and FWHM) of Ha and He I 
for the five stars are given in Table 6. A comparison with 
the EWs of Ha predicted by Kurucz models shows a general 
agreement with the UV spectral classification of Table 2. In 
the spectrum of HD 53623 both lines show a narrow, deep 
and symmetrical central absorption core. Since this object 
is a single star, it suggests that this central component is 



formed in a shell. The FWHM of the shell component of Ha 
is ~ 40 km s~^ and that of the He: line is ~ 60 km s~^. 

Apart from this shell component in the Ha and He I 
profiles of HD 53623 there is a similarity between the He I 
profiles of HD 53623, HD 53755 and HD 52721 and also be- 
tween the Ha profiles of the first two stars. This is especially 
expressed by the similarity in their half-widths and proba- 
bly also by their values of vsini. If we lay the corresponding 
profiles on top of each other we observe differences in central 
depth. Since the spectral types of these stars are comparable 
it is conceivable that (similar to HD 52721) the Ha and He: 
profiles of HD 53623 have emission components with respect 
to those of HD 53755. There is also some similarity between 
the line profiles of HD 53974 (FN CMa) and HD 54439, but 
the centre of the He I line of HD 53974 is deeper than that of 
HD 54439, which could be due to a difference in luminosity 
type or in the contribution of an emission component for 
HD 54439. We shall return to this point in the discussion 
(Sect. 5). 

The Ha profile of HD 52942 (FZ CMa) is very broad 
and has a core which can either be double or filled in by 
emission. The double nature of the profile seems most prob- 
able, since the object is a double-line eclipsing binary with a 
period of 1.27 days (Moffat & Vogt 1974) with a maximum 
radial velocity separation of 342 km s~^ (at phase 0.25 and 
0.75). From the ephemeris of Moffat & Vogt we conclude 
that our spectrum was taken at phase 0.635 from the pri- 
mary minimum. The separation of the Ha components in 
our spectrum is ~ 200 km s~^, which is compatible with 
this phase. 

Most of the UV lines are formed in the photospheres 
and cool circumstellar regions of these stars. The profiles of 
the Si IV and C iv resonance lines are asymmetrical and their 
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Figure 8. High resolution lUE profiles of the Si IV resonance lines 
for HD 52721, HD 53367, HD 53755 and HD 53974. The arrows 
near the stellar names point to the corresponding flux scales. 

blue wings extend to velocities of ~ 1600 km s~^ (Fig. 8; for 
HD 53974 see also Friend, 1990). This gives evidence for 
fast outflowing winds from these stars. Accretion discs can 
not accelerate the winds to such high velocities so that there 
must be a different driving mechanism for the winds of these 
stars. 



4 VARIABILITY 

We shall discuss here the photometric variability of the four 
Be stars and some first attempts to interpret the variations. 

4.1 HD 52721 (GU CMa) 

This star is a visual binary (ADS 5713) with a separation of 
0.65" and Am = 0.95 mag. The Hipparcos photometry (ESA 
1997) is consistent with a period of 1.61037 days. During the 
period 1987-1998 extensive photometry (320 observations 
over 12 seasons) of this star was collected at Mt. Maidanak 
(Uzbekistan) in the framework of the ROTOR photomet- 
ric program. The light-curve over 10 years does not show 
clear (> 0.1 mag) indications of long-term variations of up- 
per and lower limits in V . The results of a period search 
on these data has recently been published (Ezhkova 1999). 
Two periods with equal probability were found: 1.610158 
and 0.80508 days. For the shortest period the date of the 



minimum is given by: JDmin = 2447018.115 + 0.80508 x E 
days. Earlier monitoring of HD 52721 in UBVRI during 5 
nights (107 observations) in 1985 (Praderie et al. 1991) ap- 
pear to be in excellent agreement with these results (Fig. 10) . 
From the photometric data of 1985 we derived the extinc- 
tion corrected UBVRI fluxes at maximum and minimum 
brightness and compared the variations with the correspond- 
ing intrinsic fluxes of the Kurucz model. It appears that the 
flux variations from maximum to minimum brightness cor- 
respond roughly with a change in T^s from 25,000 K (Bl V) 
to 22,500 K (B1.5 V). The simultaneously obtained high res- 
olution Ha profiles (at phases a, b, c, d and e in Fig. 10) 
show no difference in shape but a 40% increase in EW at 
photospheric minimum c (Praderie et al. 1991). Several in- 
terpretations of this periodicity may be considered: 

a. An eclipsing binary system. The hght curve for the 
longest period is consistent with a binary system with two 
components with equal radii and effective temperatures. 
However, the radial velocities of the He I 5875A line in a 
sequence of 0.7 A resolution spectra of HD 52721, taken in 
1991 over a period of several weeks with the main spectro- 
graph on the 6-m SAO telescope, show variations with an 
amplitude not larger than « 20 km s~^, which is 10-15 times 
smaller than expected for such a system with components 
of equal mass (5-10 Mq). The radial velocity variations in 
the profiles of He I 6678 A, observed by Corporon & La- 
grange (1999) are consistent with the SAO data. However, 
if the mass of the secondary is at least five times smaller 
than that of the primary component, an upper limit of « 30 
km s~^ is not in conflict with the shorter period of 0.80508 
days. In this case we do not observe a secondary minimum in 
the visual light curve because the eclipse is only partial. Fur- 
ther support for the hypothesis that we are dealing with a 
secondary that is much less massive than the primary comes 
from the observation that the colours of the system become 
slightly redder at minimum brightness (Fig. 10) . With a dis- 
tance of 1050 pc and T^s = 25,000 K we derive a luminos- 
ity of 1.8 X 10'' Lq and a radius of 7.1 R© for HD 52721. 
From theoretical evolutionary models (Schaller et al. 1992) 
we then estimate that the mass of the primary should be 
close to 12 Mq. The orbital distance is then « 9 Rq, which 
is large enough to accommodate a 1 M© secondary star. We 
estimate that the secondary could be a KO star (if its age is 
not much larger then that of the primary). 

b. A pulsation instability. In the HR diagram HD 52721 
is situated close to the predicted border of the 13 Cep insta- 
bility zone for masses 10-12 Mq (Pamyatnykh 1999). The 
pulsation periods predicted by the model depend on vari- 
ous input surface parameters such as hydrogen and metal 
abundances and rotation. For a mass of 12 Mq and logT = 
4.35 (B1.5 V), solar abundances and no rotation, pulsations 
with periods between « 0.4 and 1.6-2.0 days are possible in 
the second order mode {I = 2). For logT = 4.3 (B2 V), the 
period can range from 0.8 up to 2.4 days. Both measured 
periods are in this range. It is not clear if the magnitude 
of the photometric and line proflle (Ha and He l) variations 
can be explained by these pulsation models. 

4.2 HD 53367 (V750 Mon) 

HD 53367 is a visual binary with a separation of 0'.'65, a 
magnitude difference of 1.41 mag and a position angle of 
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Figure 9. High resolution Ha (left) and He I 5876 A (right) profiles of HD 53623, HD 53755, HD 53974 and HD 54439. 
Table 6. Ha and He I Equivalent Width (EW) and Full Width Half Max. (FWHM) for five absorption-line stars in CMa Rl. 



Line 


HD 53755 


HD 53974 


HD 54439 


HD 53623 


HD 52942 


Ha 












EW [A] 


3.02 ± O.fl 


2.54 ± 0.18 


3.56 ± 0.10 


4.75 ± 0.06 


3.71 ± O.V. 


EW* [A] 


3.2f ± 0.09 




3.70 ± 0.19 






EW** [A] 








4.34 ± 0.12 




FWHM [km s-i] 


480 


280 


320 


460 


807 


FWHM* [km s-i] 


420 




285 






Teff-mod. 


25,000 


30,000 


25,000 


25,000 


20,000 


log 3-mod. 


3.9 


4.0 


4.14 


4.43 


3.59 


Sp. class 


Bf IV 


BO IV 


BlIV/V 


Bl V 


B2IV 


He: (5875.6 A) 












EW [A] 


1.03 ± 0.10 


1.06 ± 0.07 


0.76 ± 0.05 


0.84 ± 0.05 




EW* [A] 


1.22 ± 0.07 










EW** [A] 








0.37 ± 0.06 




FWHM [km s-i] 


400 


200 


210 


250: 




FWHM* [km s-i] 


435 











EW* = EW after estimated correction for emission component. EW** = EW without "shell" contribution (for HD 53623). For Ha the 
corrected EWs are compared with the predictions of the models of Kurucz (1991) which lead to interpolated values of logg and 
therefore serve to check the luminosity classes of Table 2. 
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Figure 10. Light curve (left) and phase diagram (right) of HD 52721 in BVRI during 5 nights in March 1985. For the construction 
of the phase diagram a period of 0.80508 days was assumed. Vertical marks indicate the times at which Ha (top) and Mgll 2800 A 
(bottom) data were obtained. 



298° (ESA 1997). HD 53367 is also a spectroscopic binary 
(Finkenzeller & Mundt 1984). The main component shows 
radial velocity variations in the He: 4471, 5876 and 6678 A 
lines, which can be fitted with a period of 166 days and an 
amplitude of about 20 km s~^ around an average velocity of 
+48.2 km s~^. A possible orbital solution has e = 0.18 and 
u) = 303° (Corporon & Lagrange 1999). 

HD 53367 shows a low vsini (20-40 km s"^). There 
are several indications that this is not due to a small incli- 
nation, i.e. that we are not observing the star pole-on. One 
indication is the observation by Oudmaijer & Drew (1999) 
of the variation of the linear polarization percentage over 
the Ha profile of HD 53367, which suggests that the ionized 
envelope, projected on the sky, is not observed as a circu- 
lar disc as expected in case the star is observed pole-on. 
Another indication against a pole-on view of HD 53367 has 
been given by Oegerle et al. (1983). These authors argued 
that the simultaneous occurrence of emission in the 8446 A 
O I line and absorption in the 1302 A O I line can only be 
understood if the star has an envelope, which is observed at 
large inclination with the polar axis (see also Sect. 3.3.1). 

Apart from its velocity variations, HD 53367 also shows 
photometric variability. Hipparcos observed an amplitude of 
0.24 mag, but no period is given (ESA 1997). However, over a 
longer interval than these Hipparcos observations, the pho- 
tometry of the ROTOR group, the Wesleyan Observatory 
and the Corralitos Observatory suggests a mean periodic 
variation with a period of about 9 years with an amplitude 
of about 0.3 mag (Fig. 11). The colour-magnitude diagram 
of HD 53367 shows that the colours are bluer at minimum 
brightness. Nevertheless after correction of the UBVR fluxes 
for extinction it appears that in the minimum phase the pho- 



tospheric temperature has decreased from 30,000 K to about 
25,000 K. 

At minimum brightness (Oct. 15, 1989) and at max- 
imum brightness (March 2, 1991) high resolution spectra 
have been taken with the Main Stellar Spectrograph on the 
6-m BTA telescope of the SAO (Fig. 7). They show that 
in the minimum phase the emission component in all the 
lines is weaker than in the maximum brightness phase. The 
emission of H/3 at minimum brightness is very faint. This 
could be related with variations in the Ha profile, reported 
by Halbedel (1989), but the time-scale of these profile varia- 
tions is not known yet, except for small changes in Ha after 
1 month (Fig. 4a). Ha profiles of May-Dec. 1974 (Garrison 
& Anderson 1977), July 1981 (Finkenzeller & Mundt 1984), 
1983 (Andrillat 1983), June 1986 (Hamann & Persson 1992), 
Nov. 1987 (Halbedel 1989) and our observations of Dec. 1994 
and Jan. 1995 (Fig. 4a) show differences, although they are 
small. However, the EWs of Ha show a strong decrease be- 
tween Dec. 1987 and Jan. 1992 (Table 4), a period for which 
no Ha line profile information is available. For H5 (Fig. 5), 
the contribution of emission to the profile is small. We find 
an EW of 3.31 ± 0.08 A for B.5 in 1991 and 3.73 ± 0.16 A 
for this line in 1989. These values are consistent with the 
predictions of Kurucz models for TcS = 30,000 K, logg = 
4.25 (BO IV) or 25,000 K and logg = 3.75 (BOIII/IV). 

It may be noted here that HD 53367 is not the only 
young Be star with a long photometric period. Another 
early type Herbig Be star, BD-|-65°1637 shows a similar 
behaviour, with a period of at least 15 years. It is to be 
expected that more long periodicity young variables will be 
detected when the photometric observations are prolonged. 
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Figure 11. Light curve (top) and colour-magnitude (bottom) 
diagrams of HD 53367. 
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Figure 12. Light curve (top) and colour-magnitude (bottom) 
diagrams of LkHo 218. 



We can think of several possible causes for such long photo- 
metric periods: 

a. Magnetic activtty of a B star. Little is known about 
the magnetic fields of early type B-stars. The intrinsic con- 
tinuum polarization HD 53367, observed by Oudmaijer & 
Drew (1999), is very low and has an angle of 44 degrees, 
which agrees with that of the dust arc in CMa Rl (Vrba, 
Baierlein & Herbst 1987). This suggests that the rotation 
axis of HD 53367 "remembers" the larger scale circumstel- 
lar field direction, which means that at least at the time 
of formation of the dust arc, the star had a magnetic field 
that could be coupled with the large scale field of CMa Rl. 
The exceptionally low rotation velocity of HD 53367 may 
therefore be a result of magnetic braking during outfiow or 
accretion. 

Attempts to measure the magnetic field of HD 53367 
have been made recently by Glagolevskij & Chountonov 
(1998). Although no significant field was detected at that 
moment it does not exclude the possibility of the presence 
of weak surface fields at other moments. We suggest that 
every 9 years during a short period the internal field moves 
upward and emerges in the surface layers of the star, thereby 
producing spot formation, similar to the solar activity cycle. 
If the spots cover a large part of the surface they may lower 
the effective temperature during the activity phase and may 
have an effect on the photometric colours. If the surface 
magnetic field was stronger during the earlier evolutionary 
phases of the star (e.g. during the convective phase) the 
field could also have braked the rotation of the star during 
outflow or accretion. 

The braking time due to angular momentum loss in out- 
fiow has been calculated by Friend & MacGregor (1984). If 
the Alfven radius Ra of the magnetosphere is larger than 
the corotation radius Rc the matter trapped in the magne- 
tosphere at Ra will lose angular momentum and will brake 
the rotation of the star. For a 20 Mq star with a rota- 
tion velocity of 250 km s~^ and a 200 G surface field we 
find an Alfven radius of 1.74 R*. The outflow velocity at 
Ra is then ~ 700 km s~^. From equation 38 of Friend & 
MacGregor's paper it is then seen that a mass outflow of 
4 X 10"'' Mq yr"^ in a field of only 200 G is able to brake 
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Figure 13. Light curve (top) and colour-magnitude (bottom) 
diagrams of LkHa 220. 



the rotation of the star within 0.6 Myr. Similar results are 
obtained if one estimates the braking during the pre-main 
sequence accretion phase (Rathnasree 1994). At an accre- 
tion rate of 4 x 10"'^ M© yr"\ a field of 400 G will be large 
enough to brake the rotation from close to break-up velocity 
to very slow rotation in about 1.7 Myr. In order to arrive at 
the same rotation velocity after a time-scale of 6 Myr smaller 
accretion rates and smaller flelds would be sufficient. 

b. A rotationally distorted Be star. The long-period pho- 
tometric variation may be due to a quadrupole contribution 
in the external potential of a rotationally distorted Be star 
(Papaloizou, Savonije & Henrichs 1992). We should then ex- 
pect a simultaneous periodic variation in the V/R ratio of 
the Balmer emission profiles. We have noted a change from 
V/R < 1 for the profile observed in 1974 to values of V/R 
> 1 for the profiles obtained in 1981 and 1994, but we do 
not have enough high resolution Ha observations to detect 
a periodicity. 
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4.3 LkHQ 220 (HU CMa) & LkHa 218 (HT CMa) 

In the course of the UBVR photometric monitoring program 
(ROTOR) at Mt. Maidanak, 153 observations of LkHa 218 
and 168 observations of LkHa 220 were made between 
1985 and 1995. The minimum time interval between the 
observations is 0.8 days. The variability was classified by 
Shevchenko et al. (1993) as LQ and BF respectively. Both 
stars vary irregularly with low amplitude 0.07 mag) on 
short time scales (2 days) but around 1990 the lightcurve 
of LkHa 220 also shows a broad and shallow minimum dur- 
ing about 600 days (Fig. 13). For both stars there were a 
few variations with larger amplitude (0.2-0.4 mag). These 
variations were characterised by reddening during a decreas- 
ing brightness, as expected from an increasing circumstellar 
dust extinction. The slope of these variations in the colour- 
magnitude diagram is close to 3. The presence of circum- 
stellar dust is also evident from the SEDs (Fig. 1). 



5 SUMMARY AND DISCUSSION 

In the previous sections we have studied the photometric 
properties and the high resolution profiles of several lines in 
the spectra of the early B-type stars HD 52721, HD 53367, 
LkHa 220, LkHa 218, HD 52942, HD 53623, HD 53755, 
HD 53974 and HD 54439, which are members of CMa Rl 
(paper I). The membership is confirmed by the common 
foreground excess. We have described the photometric be- 
haviour of the first four stars and we have discussed the pos- 
sible origin. The conclusions of this discussion can be sum- 
marized by: (a) HD 52721 varies with a period of 0.80508 
days and the variation is consistent with the presence of a 
partially eclipsing binary companion, (b) HD 53367 has a 
long-term photometric variation with a period of « 9 yr, 
and a shorter variation (166 days) in radial velocity. The 
second variation (detected by Corporon & Lagrange 1999) 
is probably due to the presence of a binary companion. We 
have given arguments for the interpretation of the long- 
term photometric variability in terms of a cycle of mag- 
netic activity. During its previous evolution, an outflow of 

5 X 10"'^ Mq yr~^ in the presence of a modest (200 G) mag- 
netic field would have been sufficient to brake the rotation of 
the star to its present low value of v sini. (c) The photomet- 
ric variability of LkHa 218 and LkHa 220 seems partly due 
to variability of the column density in their circumstellar 
dust envelopes. 

Three of these stars (HD 52721, HD 53367 and 
LkHa 220) have spectra, which are characterised by emis- 
sion components in Ha, H/3, Hel (5867 A), the Call (2) red 
triplet, the O I 7773 A triplet and the O I 8446 line, which 
gives arguments for the presence of a disc around these stars 
(see also Hamann & Persson 1992). Whether these discs are 
classical Be star discs or the remnants of pre-main sequence 
accretion discs remains uncertain so far. Note that stellar 
winds from early B type main sequence stars and classical 
Be stars usually reach velocities in excess of 1500 km s~^ (see 
e.g. Porter 1999), much higher than the 600 and 450 km s~^ 
as observed for HD 52721 and HD 53367 in the wings of 
the Si IV lines (Fig. 8) and the 340 and 370 km s"^ in the 
wing of the 1548 A Civ line. However, recently Drew, Proga 

6 Stone (1998) modelled the wind from accretion discs of 



young early type B stars and found that rotating discs can 
have high mass-loss rates (> 10~^ Mq yr~^) in combina- 
tion with lower velocity (a few hundred km s^^) equatorial 
disc winds. These are precisely the mass losses and wind 
velocities we measured for HD 52721 and HD 53367. 

We also discussed the Ha and 5876 A Hel profiles of 
the five stars without clear indications of emission in Ha. 
In these stars Ha and Hel have very similar profiles, which 
suggests that Ha is mainly formed in the photospheres of 
these stars. The high wind velocities measured in the wings 
of their Siiv profiles (« 1500 km s"^) of HD 53755 and 
HD 53974 are typical for normal early type stars. If the 
four Ha emission stars indeed have accretion discs, the next 
question is whether the other five young stars HD 52942 
(B2IVn), HD 53623 (BO.5-1 IV/V), HD 53755 (BOV + 
F5 HI), HD 53974 (BO Illn) and HD 54439 (Bl V) could have 
such discs. Although HD 53974, HD 53755 and HD 54439 
show weak indications of emission in their Ha and He I pro- 
files (Fig. 9), it seems doubtful if discs are present around 
these stars. 

One obvious question is then why these stars do not 
have accretion discs in spite of their ages being compara- 
ble with those of the stars with clear Ha emission (paper 
I). The conditions for the existence of discs could be infiu- 
enced by the presence of nearby binary companions. Both 
HD 52721 and HD 53367 have such companions, but of the 
other group also HD 53974 and HD 52942 have companions 
at comparable distances. We suggest that the survival of an 
accretion disc depends on infiuences from outside (winds or 
radiation from O stars or shock waves from nearby super- 
novae) which could evaporate accretion discs around young, 
early-type B stars (Yorke & Welz 1996). A natural candi- 
date for such an outside disturbing source is the supernova 
about 0.8 Myr ago, which accelerated the star formation in 
CMa Rl (Herbst & Assousa 1978; see also paper I). 

A strong indication for the position of this supernova 
is given by the centre of curvature of the ionization front 
S296, which is near a (1950) = 7''08™, S (1950) = -11.2° 
(Reynolds & Ogden 1978) and close to the original position 
of the runaway star HD 54662 (07III), deduced from its 
motion in space since the SN explosion 0.8 Myr ago. The 
other runaway candidate, HD 57682 (09 Ve; Comeron et al. 
1998), passes the centre at a larger distance (if the Hipparcos 
proper motions are assumed). It is remarkable that within 
the circle of curvature with radius 1^4 (« 25.7 pc) around 
the SN position only one star with Ha emission, HD 55135 
(B4 Ve) (near the border and a foreground star at a distance 
of « 500 pc with an E(B — V) of 0.10 mag) is present. This 
is in sharp contrast to the field at the west side of S296, 
where in a dust cloud region mapped with IRAS (see paper 
I and Luo 1991), the nine stars of our study are located. In 
the field between 6''50"-7''02™ and between 8°30'-12°00', 
bordering S296 at its west side, Tovmasyan et al. (1993) 
have found 23 stars at a distance of 1100 pc from which 
six stars have Ha emission. If we add LkHa 218, LkHa 220 
and Z CMa, we have 26 stars in this field from which nine 
stars have Ha emission. The remarkable contrast between 
the relative numbers of Ha emission line stars east and west 
of S296 suggests that the dust clouds near this ionization 
front shields the stars westward of S296 from UV irradiation 
by the SN or by the two 07 stars in the circle (HD 54662 
and HD 53975). 
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Recently several investigations have been devoted to the 
photo-evaporation of proplyds (discs) during FUV and EUV 
irradiation by an external star (Johnstone, HoUenbach & 
Bally 1998; Richling & Yorke 1998; Storzer & HoUenbach 
1999). Until now few numerical estimates have been made 
and mostly for low mass stars. However, Richling & Yorke 
(1998) have also computed the evolution of a 1.67 Mq disc 
around a massive (8.4 Mq) star as a result of irradiation 
from outside with a flux of 10^'^ EUV photons cm~^ s^^ 
(case A, model I) during 3000 yr. The authors do not give 
estimates for the Hq emission of this disc, as they did for the 
low mass discs, but although the 1.67 Mq disc is strongly 
deformed by the EUV influx, it seems to have conserved its 
inner structure after the period of 3000 yrs. In the CMa 
Rl situation the O stars have an EUV photon flux of 1.5 x 
10''^ s~^ (Cruz-Gonzalez et al. 1974) at most, which for a 
distance of 18 pc to the dust barrier implies an upper limit of 
3 X 10* EUV photons cm~^ s~^ at the location of the barrier 
during perhaps 6 Myr. In comparison with the numerical 
simulation of Richling & Yorke these fluxes seem to be too 
low to affect the discs of the stars east of the dust shield. 
The presence of several Hq emission stars close to HD 54662 
but outside the circle (HD 55439, HD 54858, MWC 551) 
seems to confirm this. The situation is quite different for the 
supernova. From a typical SN type II such as SN 1979c in 
MlOO (Panagia et al. 1980) a flux of « 300 ergs or 3 x 10^^ 
FUV photons cm^ s~^ was flowing through a sphere at a 
distance of 25.7 pc at the seventh day after the explosion. 
Since the decay of the luminosity is exponential, the initial 
flux may easily have been higher by one order of magnitude. 
For other type II supernovae such as the Crab SN and SN 
1181 (Panagia & Weiler 1980), the estimated FUV fluxes 
at 25.6 pc are of the same order. Such fluxes may initially 
be 300-1000 times larger than assumed in the simulation of 
Richling & Yorke, but after 200 days may have decreased 
to 3 X 10^" FUV photons cm~^s~'^ (if the decay is similar 
to that of SN 1987A (CassateUa et al. 1987), which is not a 
typical SN type II). Although the disc simulations for UV 
fluxes with such time-behaviour have not been made yet, we 
expect an increased erosion rate of the outer disc (emission) 
layers of stars in the neighbourhood of the ionization front. 
However, discs of stars westward of this front are shielded 
by extinction of the dust cloud. 

Although these cloud extinctions in the direction per- 
pendicular to the line of sight are unknown, we can make 
an estimate of these from the 100 jxm optical depth distri- 
bution in the dust cloud, derived by Luo (1991) from the 
IRAS fluxes. We assume that the stars are in the mid-plane 
of the dust cloud at a geometrical depth of 15 pc (the ex- 
tent of the cloud on the plane of the sky is « 30 pc) and 
extrapolate the optical depth to the light path from the SN 
to the star. To this cloud extinction in the light path we 
add the circumstellar extinction, determined in the line of 
sight to the star. This total absorption in the light path from 
the SN to each star gives the following percentages of trans- 
mission of the FUV (at ^ 1500 A) flux: 12% (HD 52721), 
4.5% (HD 52942), 1% (HD 53367), 1.2% (LkHa 218), 
6% (LkHa 220), 60% (HD 53623), 40% (HD 53755), 26% 
(HD 53974), 28% (HD 53456), 40% (HD 53691), < 10"-*% 
(Z CMa). 

This result shows that the stars with emission discs have 
received the lowest FUV photon fluxes from the SN, whereas 



the stars with high FUV transmissions show no evidence for 
the presence of such discs. The fact that HD 52942 shows 
no evidence for an emission disc in spite of its low incident 
FUV flux may be related to the binary nature of this ob- 
ject. We conclude that the presence of emission discs in and 
around CMa Rl appears to be correlated with the degree of 
shielding of the discs against photons from a central super- 
nova source. A full explanation of emission discs in the CMa 
OBI association will require numerical simulations of disc 
evaporation by the FUV photon flux from the supernova. 
However, other effects of the SN explosion (EUV irradia- 
tion, shock waves) are largely unknown, and should also be 
investigated. 
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